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Abstract. We prove that bounded real analytic domains in K" with the symmetries of 
an ellipsoid, and with one axis length hxed, arc determined by their Dirichlet or Neumann 
eigenvalues among other bounded real analytic domains with the same symmetries and axis 
length. Some non-degeneracy conditions are also imposed on the class of domains. It follows 
that bounded, convex analytic domains are determined by their spectra among other such 
domains. This seems to be the first positive result on the well-known Kac problem, can one 
hear the shape of a drum?, in higher dimensions. 
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1. Introduction and the Statement of results 



The purpose of this article is to prove that bounded analytic domains fl C M. n with ± 
reflection symmetries across all coordinate axes, and with one axis height fixed (and also 
satisfying some generic non- degeneracy conditions) are spectrally determined among other 
such domains. This inverse result (Theorem 1) gives a higher dimensional analogue of the 
main result of [Z2] that "bi-axisymmetric" real analytic plane domains are spectrally deter- 
mined among other bounded analytic plane domains with the symmetry of an ellipse. To 
^ ! our knowledge, it is the first positive higher dimensional inverse spectral result for Euclidean 
domains which is not restricted to balls. Negative results (i.e. constructions of non-isometric 
isospectral pairs) are given in [U, GW, GWW] (see also [GS] for some non-Euclidean do- 
mains). Higher dimensional inverse results for semi-classical Schrodinger operators with 



O ! similar symmetries have recently been proved in [GU, H]. 

1.1. Statement of results. We consider the eigenvalue problem on the domain fl with 
^ ■ the Euclidean Laplacian and with boundary conditions B: 

a; 

B<p J (y) = 0, yedfl. 

The boundary conditions could be either Dirichlet Bip = <p\dn, or Neumann Bip = d u (p\on 
where d u is the interior unit normal. 

The (Z/2Z) n symmetries of the title are the maps 

(2) Cj . (yX\, . . . , X n ) y (^-1) • • • j %ji • • • j %n) 

and we assume that they are isometries of fl. The symmetry assumption implies that the 
intersections of the coordinate axes with fl are projections of bouncing ball orbits preserved 
by the symmetries. We recall that a bouncing ball orbit 7 is a 2-link periodic trajectory of 
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the billiard flow, i.e. a reversible periodic billiard trajectory that bounces back and forth 
along a line segment orthogonal to the boundary at both endpoints. The endpoints of the 
projection to Q of the bouncing ball orbit are fixed points of all but one of the isometries Uj\ 
the remaining one fixes the projected orbit setwise but interchanges the endpoints. We add 
the generic condition that at least one of these bouncing ball orbits is non- degenerate (see 
(3) for the conditions). We also fix the length L 1 = 2L of this bouncing ball orbit 7. 

We denote by T>l to be the class of all bounded real-analytic domains Q C W 1 satisfying 
these assumptions. Thus, V L is the class of domains such that: 

(3) 

' (i) <Jj : f2 — > Q is an isometry for all j — 1, . . . , n; 

(ii) one of the coordinate axis bouncing ball orbits, called 7, is of length 2L 

(iii) the lengths 2rL of all iterates 7 r (r = 1, 2, 3, ... ) have multiplicity one in Lsp(fl); 
< (iv) 7 is non-degenerate, i.e. 1 is not an eigenvalue of its Poincare map P 7 ; 

if 7 is elliptic and {e ±tai , ...e ±M * n ~ 1 } are the eigenvalues of P 7 , we 

further require that {cci, ...,a ra _i} are linearly independent over Q. We assume the 

same independence condition in the Hyperbolic case or mixed cases. 

Here, Lsp(Q) is the length spectrum of Q, i.e. the set of lengths of closed billiard trajec- 
tories (cf. [PS, Z3]). Multiplicity one means that there exists precisely one closed billiard 
trajectory of the given length up to time reversal. Let Specs (f2) denote the spectrum of the 
Laplacian A^ of the domain Q with boundary conditions B (Dirichlet or Neumann). 

Theorem 1. For Dirichlet (or Neumann) boundary conditions B, the map Specs '■ T^l —> ^+ 
is 1-1. 

In other words, if two bounded real analytic domains Q 1? f2 2 C M. n possessing the symme- 
tries of an ellipsoid and satisfying the non-degeneracy and length assumptions of (3) have 
the same Dirichlet (resp. Neumann) spectra, then they are isometric. To our knowledge, the 
only prior positive result on the inverse spectral problem for higher dimensional bounded 
domains is that a domain with the Dirichlet (or Neumann) spectrum of a ball must be a ball 
[KAC] . In that case the proof is based on the trace of the heat semi-group rather than the 
wave group or resolvent kernel. The heat trace for the Dirichlet (or Neumann) Laplacian 
of a bounded domain has the singularity expansion, 

Tre tA % ~ t- n ' 2 (C n Vol n (Sl) + C' n V oln^dtyr 1 ' 2 + •••), * -> 0+ 

where C n , C' n are constants depending only on the dimension. Hence the volume and surface 
measure are spectral invariants. The ball is determined as the unique domain where the 
isometric inequality Vol n -i(d£l) > A n Vol n (Q)~^~ (for a certain constant A n ) is an equality. 

Our proof of Theorem 1 has a similar form in that we calculate some special spectral 
invariants and then use the invariants to uniquely determine the domain. But instead of the 

heat semi-group we use the wave group e rf "v^n or more precisely the semi-classical resolvent 
R^{k) = — (Aq + k 2 )^ 1 for k G C, which is a semi-classical Laplace transform of the wave 
group (see §2.3). Here, we are assuming that the boundary conditions are Dirichlet, but the 
methods and results are valid for the Neumann Laplacian A^ with only minor modifications. 
The spectral invariants we study are the 'wave invariants' associated to one of the bouncing 
ball orbits defined by the coordinate axes. The key advantage of these wave invariants is 
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that they are localized at the endpoints of the projected orbit, whereas heat invariants are 
integrals of curvature invariants over Q or dQ. In Theorem 2, the wave invariants of bouncing 
ball orbits are expressed in terms of the Taylor coefficients of the defining function of Q near 
the endpoints. Under our symmetry assumptions, the Taylor coefficients are determined 
from the wave invariants. That proves Theorem 1. 

As a corollary, we obtain a result for convex analytic domains that does not require any 
length to be marked. 

Corollary 1. Let C be the class of analytic convex domains with (Z/2Z) n symmetry, such 
that the shortest closed billiard trajectory 70 is non- degenerate and satisfies the conditions 
{iii) and (iv) of (3). Then Specs: C — > IR+ is 1-1. 

This follows from Theorem 1 and a result of M. Ghomi [Gh] that the shortest closed 
trajectory of a centrally-symmetric convex domain is automatically a bouncing ball orbit. 
Hence the length of this orbit is self-marked, and it is not necessary to mark the length 
L 1 = 2L of an invariant bouncing ball orbit 7. 

1.2. Balian-Bloch and wave invariants at a bouncing ball orbit. As mentioned above, 
the proof of Theorem 1 is based the study of spectral invariants of the Dirichlet or Neumann 
Laplacian of f2 known as the Balin-Bloch (or wave trace) invariants at the closed billiard 
trajectories 7 of f2. The Balian-Bloch invariants B^j are coefficients of the regularized trace 
expansion 

00 

(4) TrR° p (k) ~ V D „(k) ]T B^k-t, -+ oo, 

of the smoothed semi-classical resolvent p (k) where p is localized at the length of the 
closed orbit 7. The smoothed semi-classical resolvent is defined in (14) and the precise 
statement of (4) is given in Theorem 2.1. The factor V Dn (k) is a well-known symplectic 
factor that is reviewed in §2.3 and discussed in more detail in [GM, PS]. 

The Balian-Bloch invariants are named after the physicists who introduced them in [BB1, 
BB2] and studied them on a somewhat non-rigorous formal level. Since then, a long stream 
of mathematical works have been produced on the dual singularity expansion of the trace 
Tr costyS—Aft of the wave group. The classical results on wave trace invariants on com- 
pact Riemannian manifolds without boundary are due to Colin de Verdiere, Chazarain and 
Duistermaat-Guillemin. The wave trace expansion was then generalized to manifolds with 
boundary by Guillemin-Melrose in [GM] (see also [PS] for a very thorough study). The 
semi-classical resolvent and wave group are related by a Laplace transform (see (13)) and so 
the semi-classical (i.e. large k) expansion (4) is essentially the same as the singularity expan- 
sion of Tr cost a/— A®. Algorithms for calculating the coefficients in the boundaryless case 
were given in [G, G2, Zl]. In [Z2, Z3] an algorithm was given for calculating the invariants 
in the boundary case but it was only implemented for plane domains. In this article, the 
algorithm is developed for higher dimensional domains and explicit formulae for the Balian- 
Bloch or wave trace invariants are given in Theorem 2. The calculations also draw on the 
analysis in [H] of similar invariants for semi-classical Schrodinger operators. This result is of 
independent interest and is valid without any symmetry assumptions. 
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We now state the formulae for the Balian-Bloch invariants. They require some more 
notation which will be further discussed in §2 and §7. Almost the same notation is used in 
[Z3]. We align the axes so that the bouncing ball orbit 7 is a vertical segment of length L 
with endpoints at A — (0, ^) and B = (0, — |f), where denotes the origin in the orthogonal 
x' = (x 1 , . . . ,x n ^ 1 ) E R™ -1 plane. In a metric tube T € (AB) of radius e around 7, we may 
locally express <9f2 = <9f2 + U d£l~ as the union of two graphs over a ball B e (0) around in 
the x'-hyperplane, namely 

(5) dQ + = {x n = f + (x'), \x'\ < e}, dQ- = {x n = f_(x'), \x'\ < e}. 

We will use the standard shorthand notations for multi-indices, i.e. 7 = (71, ...7 n _i), 
ItI — 7i + ••• + In, = X^ 1 ...X^i , and by h± 2r , we mean the (n — l)-vector 

UP<1 > _ (h ll vi h 22 >P1 7,(n-l,™-l),P<?\ 
a ±,2r — V L ±,2r > a ±,2r > rt ±,2r J> 

where [h l ±^\i<i,j< n -i,i<p,q<2r is the inverse Hessian matrix of the length functional £±(x' 1 , . . . , x 
given in (11). 

The following generalizes Theorem 5.1 in [Z3] from two to higher dimensions. 

Theorem 2. Let Q be a smooth domain with a bouncing ball orbit 7 of length 2L and let 
B y r j be the wave invariants associated to Y (see cf. 2.2). Then for each r = 1, 2, . . . , and 
j there exists a polynomial P r j such that: 

(1) B rJ = P rd ({D l 7 l f + (0)},{D l 7 l f40)}) with | 7 | < 2,7 + 2; i.e. the highest order of 
derivatives appearing in B^rj is 2j + 2. 

(2) In the polynomial expansion of B^.j the Taylor coefficients of order 2j + 2 appear in 
the form {D% +2 f ± (0)}. 

(3) By,o is only a function of r , L and n, and for j > 1 

Br J = E| 7H+ i f, {(/4S)^ 2 T 2 / + (0) - (^J)^ 2 T 2 /-(0)} 

+R2rAJ 2j+1 u(Q),J 2j+1 f-m, 

where the remainder R 2r ,j(J^ 2: ' +1 /+(0), J 2 ^ +l /-(0)) is a polynomial in the designated 
jet of f ± . 

(4) In the (Z/2Z)- symmetric case, where f + = f = — /_, we have the simplified formula 

Br, = |^E| 7l=J+1 fi (^cotf yD^ 2 f(0) + R2rAJ 2j+1 W))- 

(5) In i/ie (Z/2Z) n - symmetric case, formula (4) holds with remainder in RarjiJ' 2 * f{®))- 
In the above notation, the (Z/2 Z) n -symmetry assumptions in (5) are that 



(6) /+(*') = -/-(A /±KV)) = /±0A 

where <jj denotes the reflections in the coordinate hyperplanes of IR n_1 . The first assumption 
implies that there exists a function f(x') so that the top of the domain is defined by x n = 
f(x') and the bottom is defined by x n = —f(x'). The further symmetry assumptions then 
say that / is an even function in every variable x l , 1 < i < n — 1, i.e. 
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(7) f(x\ x^ 1 ) = F^x 1 ) 2 , (a:™" 1 ) 2 ), F e C"(R n - r ). 

The value of these explicit formulae is demonstrated by applications such as Theorem 1. 
Theorem 2 may be viewed as an alternative to the use of Birkhoff normal forms methods 
for calculating wave trace coefficients as in [G, G2, Zl, Z2, ISZ, SZ]. Further discussion and 
comparison of methods is given at the end of the introduction. 

Our proof of Theorem 2 is rather different from that in [Z3, Z4]. It is based on the 
construction and analysis of a microlocal monodromy operator associated to 7, inspired by 
the works of Sjostrand-Zworski [SZ] and Cardoso-Popov [CP] (see also [ISZ]), but employing 
a layer potential analysis more closely related to that in [Z3, HZ]. The trace asymptotics are 
eventually reduced to those of a boundary integral operator in Proposition 5.6 and Corollary 
5.7, and then to the stationary phase asymptotics of a certain oscillatory integral in Theorem 
6.3. Since the method and results give a higher dimensional generalization of the analogous 
results of [Z3, Z4], there is some overlap in the arguments from the two-dimensional case; 
we have tried to minimize the overlap, but it is necessary to give complete details on the 
formulae in n dimensions since they differ in numerous ways from the two-dimensional case. 

1.3. Determining Taylor coefficients from wave invariants. To prove Theorem 1, it 
is only necessary to determine the Taylor coefficients of the defining function f — f + — — /_ 
of Q at the endpoints of a symmetric bouncing ball orbit 7 from the Balian-Bloch invariants 
given in Theorem 2 for iterates of this orbit. This is done is §7.6. The proof builds on the 
methods of [Z3, H]. 

In fact, our method could be extended to show that analytic domains with fewer symme- 
tries are spectrally determined as in [Z3, H], but for the sake of brevity we do not prove that 
here. 

1.4. Discussion and Comparison of Methods. We obtain the formulae for the wave 
invariants by applying the stationary phase method to the trace of a well-constructed 
parametrix for the monodromy operator. A secondary purpose of this article is to con- 
nect the very conceptual but somewhat abstract monodromy method of [ISZ, SZ] with the 
methods of [Zl, Z2, Z3]. The articles [Zl, Z2] implicitly used the monodromy approach 
in the form given in [BBa, L, LT]. In this article, we construct the monodromy operator 
explicitly in terms of layer potentials, using in part the methods of [CP] and in part those 
of [Z3, HZ]. The monodromy approach connects nicely with the 'Balian-Bloch' approach of 
[Z3] and simplifies remainder estimates for the Balian-Bloch (i.e. Neumann) expansion of 
the resolvent. 

In calculating the trace asymptotics, we do not put the monodromy operator into normal 
form, but rather apply a direct stationary phase analysis to the parametrix. Terms of the 
stationary phase expansion correspond to Feynman diagrams and the main idea (as in [Z3]) is 
to isolate the diagrams which are necessary and sufficient to determine the Taylor coefficients 
of the boundary defining function at the endpoints of 7 from the wave trace invariants of 
iterates of 7. In this (Z/2Z) n symmetric case, there is a unique such diagram and that is 
why the symmetries simplify the problem. It is an interesting but difficult problem to 'invert' 
the spectrum when some or all of the symmetries are absent. 
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An alternative to the approach of this article is to use quantum Birkhoff normal forms 
around the bouncing ball orbit as in [G, G2, Zl, Z2, ISZ, SZ]. It suffices to prove the 
abstract result that the quantum normal form of the Laplacian or wave group at the invariant 
bouncing ball orbit, hence that the classical Birkhoff normal form of the Poincare map, is 
a spectral invariant. At that point, one could generalize the result of Y.Colin de Verdiere 
[CV] that the classical normal form determines the Taylor coefficients of / at the endpoints 
of the bouncing ball orbit when / has the (Z/2Z) n symmetries. We plan to carry out the 
details in a follow-up to this article. The methods of this article go further, since Theorem 
2 determines much more than the classical Birkhoff normal form. 

It would be interesting to obtain a more direct connection between the the 'normal forms' 
approach and the 'parametrix approach'. In general terms, normal forms for Hamiltoni- 
ans and for canonical transformations belong to the canonical Hamiltonian formulation of 
quantum mechanics, while parametrix constructions, stationary phase methods and Feyn- 
man diagrams belong to the Lagrangian or path integral approach. Normal forms are of 
course canonical, while parametrices are not: there are many possible parametrices (finite 
dimensional approximations to path integrals), and in the inverse problem it is essential 
to construct computable ones. The two approaches are dual, and although they contain 
the same information, it is formatted in different ways. In particular, the two approaches 
highlight different features of the geometry and dynamics. 

At the present time, explicit calculations and spectral inversion for boundary problems 
have only been carried out in the Lagrangian approach, despite the existence of a quan- 
tum normal form along bouncing ball orbits [Z2]. In the simpler setting of semi-classical 
Schrodinger operators at equilibrium points, one may compare the normal forms approach 
of [GU, CVG] to the Lagrangian approach of [H]. In this inverse problem, one has a one- 
parameter family of isospectral operators depending on a Planck's constant h, whereas in the 
boundary problem one has only one operator and spectrum to work with. The Lagrangian 
calculations in [H] reproduced the inverse results of [GU, CVG], and gave stronger ones 
where some of the symmetries were removed. It directly gives formulae for wave invariants, 
which are linear combinations of normal form invariants. 

It is interesting to observe that formula in Theorem 2 is very similar to the formula in [H] 
(Theorem 2.1) for the wave invariants at an equilibrium point for a Schrodinger operator on 
]R n with a unique equilibrium point at (x,£) = (0,0). This perhaps indicates a similarity 
between the quantum normal form of the Schrodinger operator at the equilibrium point 
and that of the Laplacian at a bouncing ball orbit. The formula is also similar to a trace 
asymptotics formula of T. Christiansen for an inverse problem for wave-guides [Chr], but 
that is less surprising. 

The methods of this paper have further applications. In a work in progress [HeZ] , we use 
the wave invariants to prove a certain spectral rigidity result for analytic deformations of an 
ellipse. 

Finally, we would like to thank the referees for their suggestions on improving the expo- 
sition. 
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2. Background 

In this section, we go over the basic set-up of the problem. It is very similar to that of 
[Z3] but requires some higher dimensional generalizations. We use the same notation as in 
[Z3] and refer there for many details. 

2.1. Billiard map. The billiard map j3 is defined on 

B*dn = {(y,rj); y G dQ, rj G T*dn, \rj\ < 1} 

as follows: given (y, 77) G B*dQ, with \rj\ < 1, let (y, () G S*Q be the unique inward-pointing 
unit covector at y which projects to {y,rj) under the map Tg n Q — > T*dQ. Then follow the 
geodesic (straight line) determined by (y, () to the first place it intersects the boundary 
again; let y' G dQ denote this first intersection. (If \r]\ = 1, then we let y' = y.) Denoting 
the inward unit normal vector at y' by u y ', we let (' = ( + 2(( ■ iy)ty be the direction of 
the geodesic after elastic reflection at y', and let rj' be the projection of (' to B*,dVl. Then 
we define 

P(y,v) = iy',v')- 

The billiard map is a symplectic, hence measure preserving, map with respect to the standard 
symplectic form on T*dVl. We denote its graph of f3 by 

(8) C7 billiard := graph (5 = {(/3(z), z) \ z G B*dQ}. 
In the case of a convex domain, 

(9) C billiard = T d := { (y, -V y d(y, y'),y', V y >d(y, y'))}, 

i.e. the Euclidean distance function d(y, y') is a generating function for /3. For non-convex 
domains, this graph is larger due to 'ghost' billiard trajectories which exit and re-enter Q 
but satisfy the reflection law of equal angles at each intersection point. Such ghost orbits 
are the price one pays for using a parametrix and distance function d(y, y') defined on the 
ambient space W l . 

2.2. Length functional. We define the length functional on (dVL) m (the Cartesian product), 
by 

(10) L(y u • • • , y m ) = 1 2/i — 2/2 1 H h \y m -i ~ Vm\ + \Um ~ 2/i | - 

Then L is a smooth away from the 'large diagonals' A pp+1 := {y p = y p+ \}. The condition 
that dL = is the classical condition that each 2-link defined by the triplet (y p -i,y p ,y p+ i) 
makes equal angles with the normal at y p . Hence a smooth critical point corresponds to a 
closed m-link billiard trajectory. See for instance §2.1 of [PS]. 

2.2.1. Length functional in graph coordinates near the iterates of a bouncing ball orbit. A 
bouncing ball orbit 7 is a 2-link periodic trajectory of the billiard flow, i.e. a reversible 
periodic billiard trajectory that bounces back and forth along a line segment orthogonal to 
the boundary at both endpoints. As in the Introduction we orient f2 so that AB lies along 
the vertical with A = (0,^),B = (0, — -|). We parameterize the boundary locally as 

two graphs x n = f±(x') over the x'-hyperplane. Thus, in a small tube T e (j) around AB, the 
boundary consists of two components, which are graphs of the form y = f+(x') near A and 
y = f-(x') near B. 
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We then define the length functional in Cartesian coordinates for the two possible orien- 
tations of the rth iterate of a bouncing ball orbit by 

2r 

(11) C±(x[, x' 2r ) = \l (4+i _ x 'p) 2 + (/™±(p+i)(4+i) _ fw±(p)( x 'p)) 2 - 

p=l 

Here, w± : Z 2r — > {±}, where w + (p) (resp. u>_(p)) alternates sign starting with w + (l) = + 
(resp. W-(l) = — ). Obviously the point (x' 1: . . . , x' 2r ) = (0, ...,0), corresponding to the 
r-the iteration of the bouncing ball orbit, is a critical point of C±. 

2.3. Resolvent and Wave group. By the Dirichlet Laplacian A® we mean the Laplacian 
A = YTj=i §^ with domain {u G Hq(H) : Am G L 2 }; thus, in our notation, A^ is a negative 

operator. The resolvent R®{k) of the Laplacian A^ on Q with Dirichlet boundary conditions 
is the family of operators on L 2 (Q) defined for fceC, > by 

R°(k) = -{A° + k 2 )-\ %k > 0. 

The resolvent kernel, which we refer to as the Dirichlet Green's function G®(k,x,y) of 
Q C R™, is by definition the solution of the boundary problem: 

f (A x + k 2 )G°(k,x,y) = -5(x-y), (x,yeQ) 

(12) { 

[ Gg(k,x,y) =0, x G dQ. 

The discussion is similar for Neumann boundary conditions except that its domain is {u G 
H 2 (Q) : d u u\on = 0}, and d u G^(k, x, y) — for x G dfl, where d v is the interior unit normal. 
As in the introduction, we superscript functions of the Laplacian with B to denote either 
boundary condition. 

The resolvent with either boundary condition may be expressed in terms of the even wave 
operator E^(t) = cos(£a/— A^) as 

• POO 

(13) R^(k) = ^J e* kt E*(t)dt, (Zk>0) 

In this paper we will consider the resolvent R^(k) along the logarithmic ray k = A + ir log A, 
where A > 1 and r G R+. 

Given p G C^°(M + ), we define the smoothed resolvent R^ p {k) by 

(14) Rl p {k) :=p*{pR»{p)) = [ p{k - p){pR*{p)) dp. 

We note that p(k — p) is well-defined since p is an entire function. Let us discuss in what 
sense the integral in (14) is defined. We notice that since p G R , we have defined the 
resolvent Rq(p) by Rq(p + i0 + ). Hence we can write 

R E, P ( k ) = hm e ^ 0+ j R p{k - p)pR^{p + ie)dp 

= - lim e ^ + J k R p(p')(k - p')R^(k + ie- p') dp' 
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where the last equality is obtained by taking an appropriate contour. Thus we can also take 
this last integral as our definition of the smoothed resolvent. 

Now by (13) we can rewrite Rq (k) in terms of the wave operator as: 

RE, P (k) = i^J R p(n)e^EE(t)dtdn 
if °°p(t)e ikt EE(t)dt 

= i(p(* + V=^+A>(* ->/=*§))■ 
The Poisson formula for manifolds with boundary [AM, GM, PS] gives the existence of a 
singularity expansion for the trace of Eq (t) near a transversal reflecting ray. If we substitute 
this singularity expansion into the trace of (15) we obtain an asymptotic expansion in inverse 
powers of k for the smoothed resolvent trace: 

Theorem 2.1. [AM, GM, PS] Assume that 7 is a non- degenerate periodic reflecting ray, 
and let p G C^°(L 7 — e, L 7 + e) 7 such that p = 1 on (L 7 — e/2, L 7 + e/2) and with no other 
lengths in its support. Then for k = A + ir log A with r G R + , TrR^ (k) admits a complete 
asymptotic expansion of the form (4) where 

• V D ^{k) is the symplectic pre-factor 

V Djl {k) = C e D ( 7 )^=====; 

v /|det(/-P 7 )| 

• P 7 is the Poincare map associated to 7; 

• £d{i) is the signed number of intersections 0/7 with dQ; 

• m.y is the Maslov index of 7; 

• Co is a universal constant (e.g. factors of2n). 

Definition 2.2. The coefficients B y j are called the wave trace invariants (or Balian-Bloch 
invariants) associated to the periodic orbit 7. 

As emphasized above, the discussion in the case of the smoothed Neumann resolvent 
Rn tP (k) is essentially the same. 

3. Resolvent and the layer potentials 

The method of layer potentials ([T] II, §7. 11) solves (12) in terms of the 'layer potentials' 
G (k, x, y), d Uy G (k, x, y) G V{VL x dfl), where v is the interior unit normal to Q, and 
d v — v • V, and where Go(k,x,y) is the 'free' Green's function of M n , i.e. of the kernel of 
the free resolvent — (A + k 2 )^ 1 of the Laplacian A on R n . A key point, first recognized in 
[BB1, BB2] and put on a rigorous mathematical basis in [Z3, HZ], is that the layer potentials 
are semi-classical (i.e. non-homogeneous) Lagrangian distributions in the k parameter, with 
additional homogeneous singularities on the diagonal. In effect we wish to make use of the 
first property and suppress the second. This will be explained in §3.1. 

First let us briefly recall the method of layer potentials. The free outgoing Green's function 
in dimension n is given by: 

(16) G (k,x,y) = ^™- 2 (2vrA;|a;-y|)-("- 2 )/ 2 Hai 1 / ) 2 _ 1 (A;|x-y|). 



(15) 
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In general, the Hankel function of index v has the integral representation ([T], Chapter 3, 
§6) 

(17) Ha«(*) = (^) 1/2 e Vi/T /o°°^ 1/2 (l ~ i- z r 1/2 ds. 

The single and double layer potentials, as operators from the boundary dfl to the interior 
Q, are then respectively defined by 

f S£(k)f(x) = J dn G (k,x,y)f(y)dS(y), xeO, 

(18) { 

[ V£(k)f(x) = J dn 2d Uy G (k,x,y)f(y)dS(y), xeSl. 

Similarly, for a function / on dfl, the boundary layer potentials S(k) and N(k), as operators 
from the boundary dQ to itself are denoted by: 

f S(k)f(x) = J m G (k, x, y)f(y)dS(y), x G Sfi, 

(19) { 

[ N(k)f(x) = J m 2d Uy G (k,x,y)f(y)dS(y), x G dQ. 

Given any function g on M n \<9fi and any x G dQ, we denote by g+(x) (resp. g~(x)) the 
limits of g(w) as to — > x G 9fi from u> G f2 (resp. u> G M n \f2). The layer potentials and the 
boundary layer potentials introduced above are related by the following 



(S£(k)f) + (x) = (S£(k)fUx) = S(k)f(x), 

(20) 

(V£(k)f)±(x) = (± I + N(k))f(x), (The jump formula). 
Using the above notation we also have the following interesting formula of Fredholm-Neumann 

(21) R°(k) - Ro(k) = -V£{k){I + N(k))~ l S£(k) t , 

where Ro(k) is the free resolvent. This is true because the kernels of both sides of the equation 
are solutions to the Helmholtz equation and also the restrictions to the boundary of these 
kernels are the same by the jump formula (20). The formula follows by the uniqueness of 
the solutions of the Helmholtz equation. 

3.1. Structure of the operator N(k). We now state the precise sense in which N(k) 
is a semi-classical Fourier integral operator quantizing the billiard map of dfl when k = 
X + ir log A and A = 3lk — > oo. It additionally has homogeneous singularities on the diagonal. 
The discussion here is adapted from [HZ]. Note that here our semiclassical parameter is | 
which is a complex parameter. 

We denote S*d£l = {(y,rj) G T*dil; \rj\ = 1} and define the diagonal set 

(22) A s * d n = {(z, z); ze S*dVl} C T*dtt x T*<9fl 

Proposition 3.1. [HZ] Assume that Q is a smooth domain. Let U be any neighborhood of 
As*an- Then there is a decomposition N(k) as 

N(k) = N (k) + Ni(k) + N 2 (k), 

where N (k) is a pseudodifferential operator of order —1 , Ni(k) is a semi- classical Fourier 
integral operator of order zero associated with the canonical relation T d (cf. (9)), and N 2 (k) 
has operator wavefront set contained in U. 
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This proposition is valid whether or not Q is convex, but in the convex case Yd = Cbmiard 
is the graph of the billiard map (cf §2.1). 

From the integral formula (17), the Hankel function Ha^ 2 _ 1 (2) is conormal at z = and 

as z — > oo, 

oo 

(23) a(z) := e~ iz Hafl^z) ~ £ a 3 z~K 

3=0 

The kernel of the single-layer potential then has the form 

S(k,x,y) = G (k,x,y) = Ck n ~\k\x - y\)-W 2 Ha^^z - y\). 
For the double layer potential it follows from the identity 

^^\z) = -^\z)-^l 1 {z), 

az z 

that 

N(k, x,y) = 2 d Uy G (k, x, y) = Ck n ~\k\x - y\)-^' 2)/2 Ha«(fc|x - y\) (^-, v y ). 

\x y\ 

We now introduce cutoff functions as in [HZ]. When k = A + ir log A, we put 

1 = v?i(k - y\) + <PoW*\x - y\) + <p 2 (\x - y\, \k\) 

where tpi(t) is supported in t > t for some sufficiently small t > (see [HZ] for the choice 
of this t )> an d <Po(t) is equal to 1 for t < 1 and equal to for t > 2. (The power 3/4 in ip 
could be replaced by any other power strictly between 1/2 and 1). We define the operator 
Ni(k) to be the one with the kernel (pi(\x — y\)N(k,x,y) which has the form 

(24) N,(k, x, y) = Ck^e^-y^x - y\) ai {k\x - y\) (^^, i/ y >, 

\x y\ 

where a\{z) = z~^~a(z) has an expansion in inverse powers of z as z — > oo (see (23)), with 
leading term ,2~( n ~ 1 )/ 2 . The operator N±(k) is manifestly a semiclassical FIO of order with 
the phase function d(x,y) = \x — y\, and thus associated with the canonical relation T d (cf. 
(9)). The operators N (k) and N 2 (k) are constructed similarly from the cut offs (p and (f 2 
respectively. See [HZ] for the details. 

4. MONODROMY OPERATORS AND BOUNDARY INTEGRAL OPERATORS 

It simplifies the resolvent trace calculation considerably to reduce it to a boundary trace. 
In [Z4, Z3] this was done by taking the direct sum RR, p (k)®Rnc jP (k) of the interior Dirichlet 
and exterior Neumann resolvents, and verifying that 

Tr(R° p (k) © R& >p (k) - RoM) =P*^ lo S det(/ + N(k)). 

The construction for the interior Neumann resolvent and exterior Dirichlet resolvent is essen- 
tially the same and is omitted for brevity. This relation was stated in the physics literature, 
and we refer to [Z4] for references. In this article, we take a related but somewhat different 
method to reduce the trace to the boundary using monodromy operator ideas similar to those 
of Cardoso-Popov [CP] and of Sjostrand-Zworski [SZ]. A novel feature is that we relate the 
monodromy operators to the boundary integral operators N±(k) and N Q (k). 
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Before going into the details, let us note some motivating ideas. First, both the monodromy 
operator M(k) and the boundary integral operator Ni(k) are quantizations of the billiard 
map in the sense of being semi-classical Fourier integral operators whose canonical relation 
is Cbiiiiard- This suggests that they must be closely related. However, M(k) is microlocally 
constructed while N ± (k) is global. Further, Ni(k) is just a piece of N(k), which has the 
more complicated structure described in §3.1, and in the boundary reduction N(k) is the 
primary object. On the other hand, there is a simple exact formula for N(k) while M(k) 
is only known through microlocal conjugation to normal form. Hence our purpose here is 
to construct a monodromy operator resembling M(k) using N(k) and the layer potentials. 
In doing so, we follow the approach to monodromy operators of [CP]. Since the resulting 
monodromy operator does not seem to arise from the abstract set-up of Grushin reductions, 
the proof that the trace reduces to the boundary is not the same as in [SZ]. 

4.1. Definition of the monodromy operator. In this section we will use the basic termi- 
nologies in semi-classical analysis such as semi-classical pseudodifferential operators, Fourier 
integral operators and semi-classical wave front sets. We refer to [GS, EZ, SZ, Al] for all 
the definitions and properties. We just mention that in the following sections when we write 
T{k) ~ S(k), for two operators T(k) and S(k), we mean T(k) — S(k) is a negligible (residual) 
operator in the sense that its kernel is of order 0(k~°°) in all C s norms. 

Let 7 be an m-link periodic reflecting ray, with vertices at Vj in dQ, where there are m+1 
vertices and i> m +i — v i- Let i]j be the projection to B*.dQ of the direction of the ray 7 where 
it hits the boundary at Vj. We denote 

dl ■= {(yi,r)i),...{v m ,r) m )} 

Let Tj be microlocal neighborhoods of (vj, rjj) in B*dVt = {(y, r))\y E dQ, r\ e T*dVL, \q\ < 1}. 
We then define a microlocalization to Tj of the double layer potential operator by 



where xr 3 is a microlocal cut off supported in Tj. More precisely Xr^k) = Opi(a,j(y,r))), 

where supp(aj(y, 77)) C Tj and aj(y,r)) = 1 in a neighborhood of (vj,r)j). Here (/ + A^ (A;)) _1 
is a parametrix for I+N (k). Notice a parametrix exists because N (k) is a pseudodifferential 
of order — 1 . 

Clearly, by the jump formula (20), for any u G C(dQ), 



Pi(k) = H t {k) - H t+1 (k)r dn Hi(k) + ■■■ + (-l) m_1 iJj + ( m _i) (k)r dn Hi + ( m - 2 ) (k) ■ ■■r dn H i {k), 



where tqq : C(Q) — > C(dQ) is the operator of restriction to the boundary. Also notice by 
our notation the indices j + m and j are identified. Then we define 
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(26) p(k) = -j2m)- 

t=l 

Proposition 4.1. We have: (-A - k 2 )P(k)u = 0. 

Proof. In fact for every i we have (—A — k 2 )Pi(k)u = 0. We observe that Pi(k) = V£(k)Qi(k) 

with 

(27) 

Qi(k) = *r, - *r l+2 [(/ + N(k))* Ti+1 ] +■■■ + (-l)™" 1 ^^ [(/ + A^))*r i+m _ 2 •••(/ + W(fc))*rJ . 

The statement follows since T>£(k) maps C(<9Q) into the solutions of the Helmholtz equation. 

□ 

We now make a couple of useful technical observations: 
Proposition 4.2. We have: 

r dn H,(k) ~ (J + W (AO + Wi(A;))tfiy 

Proof. By the jump formula we know that rgnHj(k) — (I + N(k))^r r The missing term 
N 2 has its wave front set contained in U x U where U is a small neighborhood S*dtt. But 
then WF'(N 2 ) o WF'(^ T .) = 0. So this term of the composition is negligible in k. □ 

Proposition 4.3. For all 1 < j < m, we have: 

^ Vj+1 [(I + N(k))^ Tj •••(/ + N(k))V ri ] ~ * r . +1 [N^r. ■ ■■N 1 (k)V ri ]. 

Proof. We argue inductively. For j = 1, using the proof of Proposition 4.2 we have 

^r 2 [(/ + N(k))* ri ] ~ *r 2 ^i(^)*ri + ^r 2 [(/ + N (k))* ri ] ■ 

But the last term in the expression above is negligible in k, as the two semi-classical pseu- 
dodifferentials ^r 2 an d [(I + N (k))ty ri ] are microsupported in the two disjoint open sets 
T 2 and Ti respectively. Now assume the statement is true for j — 1. Then we write 

* r . +1 [(/ + N(k))V Tj •••(/ + N(k))V ri ] ~ ^r 3+1 [N^k)^ • • • N^k)*^] 

+^r 3+1 (J+iV (A;))^r 3 [N&r^ • • • Ar^rJ ■ 
Similarly, \P r (J + iV" (/c))\l/r., is negligible in fc and the second term above is negligible. □ 

Next we have the following important proposition which implicitly defines the monodromy 
operator: 

Proposition 4.4. We have 

ranP(k) ~ I + M(k), (microlocally near dj) 

where 

(28) M(k) = - X)(- 1 ) m " 1 ^i(*)*r^_ 1 ---^i(fc)*r i . 
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Definition 4.5. The monodromy operator for the m-link periodic reflecting ray 7 is the 
operator M(k) on L 2 (dQ) defined by (28). 



Proof. We define 



Mi{k) = (-ir- 1 7V 1 (fc)*r i+m _ 1 ---iViW*r i , 



- i+m— 1 

and we show that for every % 

(29) r dn Pi{k) ~ xr t + Mi(k), (microlocally near cfy). 

Since 

^ m ^ m 

— Mj(fc) = M(k), and — 7 xr • = -f (microlocally near $7), 
i=i i=i 

taking averages of equations (29) over all % implies Proposition 4.4. So we only prove (29) 
for i = 1. Since Pi(/c) = V£(k)Qi(k), by the jump formula (20) and by (27), we have 

TO— 1 

r dn Pi(k) = (/ + JV(AO) ^(-l)% r , +1 [(/ + iV(^))*r, •••(/ + N(k))* ri ] , 

3=0 

where the 0-th term of the sum is defined to be \I/rv Now we apply Proposition 4.3 to each 
term of the above sum and we get 

TO— 1 

r d nPi(k) ~ (J + N(k)) ^(-l)^r J+1 [N^k)^ ■ • • N^k)*^] ■ 

3=0 

Hence by substituting I + N = (I + N ) + JVi and collecting the terms corresponding to the 
same number of iterations of the billiard map (3 (this number is the same as the number of 
factors Ni in each term) we obtain 

TO— 1 

lanPi(k) ~ ^(-ir 1 (/-xr, +1 )[iViW*r,---iVi(A;)*r 1 ] + ( X r 1 +M 1 (k)). 
3=1 

The j-th term of the sum above is a fc-FIO corresponding to ft. We show that each of these 
terms is microlocally equivalent to near dj. Hence only xr t corresponding to ft and M\{k) 
corresponding to (3 m survive and the proposition follows. 
Let us discuss why 

(/ - X r j+l ) [Wi(fc)tf r , • • • Ar i(^)^r 1 ] ~ 0. (microlocally near #7) 

This is true because 

WF'(N 1 (k)^ Tj ---N 1 (k)^ ri ) C {(ft (y, V ), (y, V )); (y, 77) G C T j+1 x T t , 

and because / — Xr j+1 is micro-supported away from (vj,r)j). Thus by passing to a smaller 
microlocal neighborhood of in B*dQ, we obtain a negligible operator. □ 
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4.2. Microlocal parametrix for the interior Dirichlet problem in terms of M{k). 
In this section we construct a microlocal parametrix for the Dirichlet resolvent of Q near 7 
in terms of the monodromy operator M(k). It is a microlocal version of the global formula 
(21) of Fredholm-Neumann. The discussion in the Neumann case is similar. 

Proposition 4.6. Microlocally in a neighborhood ofj x 7 c T*Q x T*Q, we have 

Rg(k) - R (k) ~ -P(k)(I + M(A;))- 1 5£(A;) t . 

Proof. Let us look at the Schwartz kernels of the both hand sides in a microlocal neighbor- 
hood of 7 x 7. We show that 

Gg(k,x,y)-G (k,x,y) ~ -(P(A;)(/+M(A;))- 1 )(5£(A;)*(A;,a;,?/)), (microlocally near 7x7) 

where the operator P(k)(I + M(k))^ 1 : C(dQ) — > C(fl) acts on the first component x of the 
kernel S£(kY(k,x,y). To prove this, since microlocal solutions are unique, it is enough to 
show that for all y the right side is a solution of the Dirichlet problem 

(-A x - k 2 ) (P(k)(I + M(k))- 1 ) (Se(kY(k, x, y)) = 0, 

rdn(P(k)(I + M(k))~ 1 )(S£(k) t (k,x,y)) ~ —(Gq — G )(k, x,y), microlocally for x near d^. 

But this is clear because by Proposition 4.4, microlocally for x G dfl near d'-f, we have for 
all y G Q 

r m (P(k)(I + M(k))- 1 ) (Se(kY(k, x, y)) ~ (I + M(k))(I + M(A;))- 1 «S^(A;) t (A;, x, y) 

~-(G%-G )(k,x,y). 
We note that the operator / + M(k) is invertible by Lemma 5.5. □ 

5. Trace formula and Monodromy operators 

Here we assume L 7 is the only length in the support of p G and p(t) = 1 near L 7 . 
From Proposition 4.6, we immediately have a reduction of the trace to the boundary. 

Proposition 5.1. We have: 

Trn{Rn, P {k) ~ Ro, P (k)) ~ - Tr dn (p * k S£(kfV£(k)Q(k)(I + M{k)y v ). 

Proof. First of all we note that the regularized trace TrQn(R® p {k) — i? , p (/c)) can be microlo- 
calized to 7, i.e. if x 7 is a microlocal cutoff around 7, then 

Tr m {R^ p {k) - Ro >p (k)) ~ Tr dn (x-y(RE, P (k) - Ro, P ( k ))x,)- 

For a proof of this fact, see [Z4] §3.3. Now by Proposition 4.6, we have 

Tr Q (RR, p (k) - R , P (k)) ~ -Tr n (p * kP(k)(I + Mik^Sm 1 ) 

~ -Tr an (p* A;^(A;)'P(A;)(/ + M(A;))- 1 ). 
The formula follows by substituting P(k) = V£(k)Q(k). □ 

This formula is useful but somewhat unwieldy. As proved in [Z3], S£(kyV£(k) = Dq + Di 
where Dq is a /c-pseudodifferential operator and where D\ quantizes /3. The same analysis 
could be used here. But it is simpler to use the alternative in the next section. 
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5.1. Interior plus exterior. If we take the direct sum of the interior Dirichlet and exterior 
Neumann resolvents, then the trace formula simplifies in that we can sum up the interior 
and exterior S£(kY oT>£(k) operators to obtain -^N'(k) := ^(d/dk)N(k). 

Proposition 5.2. We have: 

Tr Rn (R° p (k) © R» >p (k) - R , P (k)) = - f p{»)Tr(N'(k - p)Q(k - p)(I + M(k - p)y l )dp. 

J R 

Proof. We first derive an analogue of Proposition 4.6 for the exterior Neumann problem. 
We then take the trace of the direct sum of the interior Dirichlet and exterior Neumann 
resolvents. 

We construct a parametrix for the exterior Neumann problem by a modification of the 
method used for the interior Dirichlet problem. The discussion of the interior was motivated 
by the double layer representation for the interior Dirichlet Green's function. For the exterior 
Neumann problem we use the single layer representation of the exterior Neumann Green's 
function, Rg c (k) - R (k) = -S£(k)(I + N\k))- l V£(k)\ where the superscript t denotes the 
transpose. This formula is proved by expressing the left side as S£(k)ip for some ip, taking 
the normal derivative from the exterior and solving for if). We then consider 

r n c(R* c (k) -Ro(k))r n o, 

where rx is the charateristic function of X. We observe that this operator is symmetric, i.e. 
equals its transpose. It follows that 

(rnc(i#c(A;) - Ro(k))rnc)(x,y) = -(r Q cV£(k)(I + iV)- 1 ^)'^) (y, x). 

Therefore, at least on the diagonal, we can use the same parametrix formula in the exterior. 

We now complete the proof of Proposition 5.2 by taking the trace on W 1 of the direct sum 
of the two operators, 
(30) 

( r n (R°(k) - R (k))r n ~ -r n P(k)(I + M{k)Y l S£{kyr n (microlocally near 7) 

\ r n c(R^ c (k) - R {k))r n c ~ -r n .P(k)(I + M(A;))- 1 5£(A;)V n c (microlocally near 7) 

In taking the trace we may cycle S£(k) t rn to the front in the first trace and S£(k) l rQc to the 
front in the second trace. We then add them to get, 

(31) Tr R n (R°(k) © R% c (k) - R (k)) ~ -Tr dn S£(kYV£(k)Q(k)(I + M(k))~ 1 , 
where S£(k) 1 o V£(k) has the kernel 

(32) / G {k,x,w)d Uy G {k,w,y)dw = —N'(k,x,y). 

For a proof of this simple fact see equation (19) of [Z5]. This indeed is why the interior 
Dirichlet and exterior Neumann problems were combined and explains the sense in which 
they are complementary. We then convolve with p. □ 

REMARK 5.3. We notice that here the exterior trace Tr (r^ c (i?^ c p {k) — i?o,p(^)) r n c ) is neg- 
ligible in k and it is added only to simplify the expression in Prop 5.1 to the more convenient 
expression in Prop 5.2. 
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We now use Proposition 5.2 to obtain asymptotics of the trace. The next step is to expand 
(7 + M) _1 in a finite geometric (Neumann) series with remainder. We have 

no 

(33) (I+M)- 1 = ^(-l) n M n + (-1)"« +1 M n " +1 {I+M)-\ 

n=0 

The following proposition shows that, in calculating a given order of Balian-Bloch invariant 
£? 7 j, we may neglect a sufficiently high remainder of the expansion (33). 

Proposition 5.4. Assume that k = A + irlog A. For each order \k\~ J in the trace expansion 
there exists n (J) such that 

(i) Tr j R p{fi)M{k - fi) n ^ +1 {I+M{k - fiy^N'ik - fi)Q{k - fi)dp = 0(\k\- J - 1 ), 

(it) Tr n RE, p (k) = Y^t\-l) n Tr J R p(/i) M(k - p)»N'(k - ix)Q{k - fidn + 0(|fc|- J - 1 ). 

Proof. Part (ii) is easily proved by combining part (i), Proposition 5.2 and (33). It remains 
to estimate the remainder and show pat (i). For this, we need to establish an L 2 - norm 
estimate for the operator norm of M(k — fi) = M(A — + ir log A) for sufficiently large r. 

The proof is implied by the following norm estimate, which is analogous to Lemma 6.2 of 
[SZ]. Let t be the constant in §3.1. We note that the monodromy operator depends on a 
choice of t although it is not indicated in the notation. 

LEMMA 5.5. Let k = A + ir log A. For every a > there exists r > and constants b, C > 
such that, 

\\M(k-fi)\\ L 2<C\k\- a <fi> b . 
To prove the Lemma, we observe that by Proposition 4.4, 

(34) \\M(k-ii)\\<C(\k\<»>) i \\N 1 (k-ri\r, 

For some integers C and b. We will not relabel these constants in the course of our estimates. 
We recall (see (24)) that Ni(k — fi) has Schwartz kernel 

(35) C{k - /.r-V^I^Wlx - y\) ai ((k -p)\x- vIXj^T, v v ), 

\x y\ 

where <fi(t) is supported in t > t for some t > 0. We estimate its norm by the Schur 

estimate, 

(36) 

WN^k < C\k- /il"" 1 sup xean J m le^+^^-y^idx - y\) ai ((k - fj)\x - y\)\dS(y) 

< C\k- ^-l e -rlo g Xt supx&m _ y | )ai(( ^ _^\ x _ y\)\dS{y) 

< C < /i > 2n e- (r *°- e)losA , 

where we estimate |A; - //| n_1 sup x ,veoci,\x-y\>to \a x {{k - n)\x - y\)\ < C\k\ 2n < fi > 2n . 

Since we can choose any small e > and also any large r > 0, it is clear that, for any 
a > 0, there exist e and r such that \\M(k — fi)\\ < C\k\~ a < /i > b . This proves the Lemma 
and hence the first part of Proposition 5.4. 

□ 
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5.2. Trace for the iterations of a bouncing ball orbit. We now analyze the trace in 
part (ii) of Prop 5.4 when it is specialized to the rth iterate y of a bouncing ball orbit, which 
has m = 2r links. We observe that Q(k — //) is a sum of terms quantizing (3°, (3 1 , . . . , f3 2r . 
Let us write q^ik — p) for the term quantizing (3f We note that based on this notation, we 
have q2r(k — /-0 = M{k — p) where M(k — p) is the monodromy operator for y. It follows 
that 

N'(k- p)Q(k- p) 

is a sum of terms quantizing f3°, . . . ,/3 2r+1 . On the other hand if we use the monodromy 
operator for y ', then (I+M(k—p))' 1 is a sum of terms quantizing /3 2r , /3 4r , • • • . Therefore 
only three terms of N' (k— p)Q(k — p)(I + M(k— /i)) -1 are associated to 7 2r and can contribute 
to the trace: 

(1) -N' Q (k - n) M(k - fx); 

(2) N>(k - v) q 2r (k -fi) = N>(k - M(k - Li); 

(3) N[{k- Li)q 2r -i{k- Li). 

Here, we use that iVp is associated to (3° and N[ is associated to (3. We notice the terms (1) 
and (2) cancel and hence only the term (3) contributes to the trace. Also we notice that for 
the r-th iteration of a bouncing ball orbit we have only two vertices and therefore for % odd 
we have Tj = r x and for i even we have Tj = T 2 . Let us denote by 

r+ = r\, and r_ = r 2 , 

the microlocal neighborhoods corresponding to the top and bottom vertex respectively. Thus, 
by these notations we have 

Proposition 5.6. Let p e C£°(R) be a cut off satisfying supp p n Lsp(Q) = {rL 7 }. Then 
(37) 

Tr(RE, p (k)) ~ -Tr J R p(p)N[(k - p)q 2r -i{k - p) dp 
~ I J2± Tr I R p(^)N[(k-p) 

(2r-l) times 

x (N^k - p)%~J(N^k - p)V r± ) ■ • • (N^k - p^jiN^k - p)^ r J dp. 

We now express this trace as an explicit oscillatory integral. We consider both principal 
(we will define the principal terms in the next section) and non-principal terms. All terms 
arise as composition of 2r Fourier integral operators quantizing /3, hence may be expressed 
as compositions of 2r oscillatory integrals. We recall that — (I + ^o) -1 Xiy Next we 
expand 

(I + N^k))- 1 = I + N^(k), 

where N-i(k) is a (— l)st order pseudo-differential operator. If we plug (/ + N-i(k))xr± f° r 
^/r ± into the expression (37), after expanding we get 
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Corollary 5.7. Let k = A + irlog A. Up to 0(\k\~°°), the trace Tr(R^ (k)) is a sum of 
2r oscillatory integrals of the form 

jo (U\ — f f f p i[n(t-C(yi,...,y2 r ))+kC{yi,...,y2r)] 
I 2r,p\*>) — Jr Jr J(dn) 2r e 

A a 2r (k- p,y 1 ,...,y 2r )p(t)dtdpdS(y 1 )--- dS (y r ) , 

where the superscript a; < a < 2r — 1, denotes the sum of the terms which contain a 
factors of N-i, and where 

C(y 1 ,...,y 2r ) = 1 2/i — 2/2 1 H h \y2r 

and <(A; -n, Vl ,..., y 2r ) e S'7 H ((^) 2r )- 

6. Principal terms 

The goal of this section is to identify the principal terms, which generate the highest 
derivative data, and to prove that non-principal terms contribute only lower order derivative 
data. 

As in [Z3], we separate out a single oscillatory integral (the principal term I 2rp ) which 
generates all terms of the wave trace (or Balian-Bloch) expansion which contain the maximal 
number of derivatives of the boundary defining function per power of k (i.e. order of wave 
invariant). 

Definition 6.1. Let 7 be a 2-link periodic orbit, and let Y be its rth iterate. The principal 
term is the term of (37) in which \l/ r± is replaced by Xr± ■ Thus, the principal term is 

(2r— 1) times 

ll, P = - J2 Tr P * N 'i( ' N ^t N iXr ± ■ ■■ N lX r T N lX r± )d/i. 
± 

This oscillatory integral corresponds to I 2rp , i.e. the one in 5.7 corresponding to a — 0. By 
Corollary 5.7, the oscillatory integral I 2 has the phase function C{y\, . . . , y 2r ) = 1 2/1 — 2/2 1 + 
' • • + \V2r — Vi\i where y p G dQ. We may write each y p in graph coordinates as (x' p , f±(x' p )). 
We will use superscripts for the n — \ components of " ). Hence the 

integral is localized to [(— e, e)" _1 ] 2r . We notice that I 2 is the sum of I 2 £~ p and I 2 ^T, where 
they correspond to the + and — term respectively. It is clear that the phase function of I 2 ^ 
is given by 

C±(x[, x' 2r ) = \l «+i - X ' P ) 2 + (^±(p+i)K+i) - f^ ± ( P )(x' p )) 2 . 
p=l 

Here, w± : Z 2r — > {±}, where w + {p) (resp. W-(p)) alternates sign starting with w + (l) = + 
(resp. W-(l) = —). 



Now we have the following Theorem 6.3. First we have a definition: 
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Definition 6.2. Let 7 be an m-link periodic reflecting ray, and let p G C^°(R) be a cut off 
satisfying supp pnLsp(Q) = {rL 7 } for some fixed r G N. Given an oscillatory integral I{k), 
we write 

TrRl p {k) = I(k) mod 0(J2^ (J 2 > /)) 

j 

*f 

TtRZJJz) - I{k) 

has a complete asymptotic expansion of the form (4), and if the coefficient of k~i depends 
on < 2j derivatives of the defining functions f at the reflection points. 

The following Theorem is the higher dimensional generalization of Theorem 4.2 of [Z3]. 

Theorem 6.3. Let k = A + irlogX. Let 7 be a primitive non-degenerate 2-link peri- 
odic reflecting ray, whose reflection points are points of non-zero curvature of dQ, and let 
p G C^°(M) be a cut off satisfying supp p fl Lsp(Q) = {rL 7 } and equals one near rL^ for 
some fixed r G N. Orient Q so that 7 is the vertical segment {x' = 0} HQ, and so that dfl 
is a union of two graphs over [— e,e] n_1 . Then 



(1) TrR^ p {k) = I^ p mod 0(£ . <(.r'/ : )) 

(2) We also have the following integral representation for I 2 in the x' p coordinates 

(38) ^2r,p = y ] / e *■ + ^ ±( - !' ' 2r ^a2r (k, x 1 , x 2 , ■ ■ ■ , x 2r )dx 1 ■ ■ ■ dx 2r , 

^ J ([-erf-* f 

where the phase £±(x[, . . . ,x' 2r ) is given in (11), and where the amplitude is given 
by: 

a 2r ' (k, x 1 , . . . , x 2r ) = C w± (x±, . . . , x 2r )A 2r ' {k, x 1 , . . . , x 2r ) + ^~Q^J^2r (^J X li ■ ■ ■ 1 X 2r) m 
Here 

A 2r ' (k, x[, . . . , x' 2r ) = rip!l 1 jai (k (x' p — x'p^i) 2 + (fw±(p)( x p ) — fw±('pu)( x 'p+i)) 2 ) 

( 39 ) 

<X p X p+1 , ^ x ' p fwJj-(p)( x p)> ifw±(p){ x p) fw ± (p+l)( x p+l)) ~t 

^( x 'p- x 'p+i) 2 +(fw ± (p)( x p)-f w± ^ + iy( x 'p +1 )) 2 J 

where ai is the Hankel amplitude in (24)- Here w + (p) = (— l) p+1 and w^(p) = 
—w + (p). Also we have identified x 2r+1 = x[. 

Proof. To prove the first part of Theorem it is enough to show that for a given a > 1, the 
coefficient of k~ j in the stationary phase expansion of I 2r , p (k), has only Taylor coefficients of 
order at most 2j — a + 1. This is shown in §5.4 of [Z3]. The second part of Theorem follows 
from the proof of Proposition 3.10 of [Z3]. It is basically just eliminating the variables t and 
p in the integral in Corollary (5.7) using the stationary phase lemma. 

□ 
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Theorem 6.3 is a crucial ingredient in the proof of Theorem 1. It gives explicit formula 
for the phase and amplitude of the principal oscillatory integrals that determine the highest 
order jet of Q in each wave invariant. The notation A^,^, a^' ± refers to the amplitude of 
the principal terms of the 2rth integral; these amplitudes contain terms of all orders in k 
and principal here does not refer to the principal symbol, i.e. the leading order term in the 
semi-classical expansion. 



7. Stationary phase calculations of IS and the wave invariants 

It is easy to see that (see Proposition 4.4 of [Z3]) we have I 2 ^ p = 7^~ and therefore 
^2r, P = 2-^3 ■ Hence it suffices to consider the + term. The oscillatory integrals I 2 * p have 
the form (38) with the phase £+ and the amplitude a^.' + . 

The only critical point occurs when x' p = for all p. We denote by Hess£ ± (0) the 
2r(n — 1) x 2r(n — 1) matrix with components 

(40) Hess£ ± (0)= | ) > hj = 1, • • • > n - 1; p, q = 1, . . . , 2r. 

It is the Hessian of C± at its critical point (x[, . . . ,x' 2r ) = in Cartesian coordinates. 

We denote by H + the inverse Hessian operator in the variables (x[, . . . , x' 2r ) at this critical 
point. That is H+ = (77ess(£ + )" 1 (0)7J, D), where D is short for (^r, • • • a^=r)- More 
precisely, 

(41) ^^EE^'fe)' 

p,q=l t,] = l V y 

Before we apply the Stationary Phase Lemma, in two subsections we state some properties 
of the inverse Hessian matrix of £+, and also some properties of the phase function C + and 
principal amplitude a!£.' + which may be derived directly from the formula in Theorem 6.3. 

7.1. Properties of Hess (£ + ) -1 . Let {vj } ±}™z\ denote the eigenvalues of the second fun- 
damental form of dfl at the endpoints of the bouncing ball orbit. Without loss of generality 
we can assume 



(42) ^,± = 7^,/ ± (0), j = l,...,n-l. 

This is because by an orthogonal change of variable (i.e. an isometry of the plane) we can 
make Hess f± a diagonal matrix. Of course when all the symmetry assumptions are satisfied, 
then Hess f± is automatically diagonal. 

The following generalizes Proposition 2.2 of [Z3]. 

Proposition 7.1. Put a jt ± — -2(l±Lv jt ±), and let A± — Diag(a j; ±) be the (n-1) x (n-1) 
diagonal matrix with the diagonal entries a,j t ± Then the Hessian 77 2r of C + at x' = has the 
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H- 



2r 





I 





I ' 


I 


A_ 


I 





< 


I 


A + I 





> / 








A- . 



where there are 2r x 2r blocks and each block is of size (n — 1) x (n — 1). 

Proof. There are 2r sets of variables rc^ and therefore there are 2r x 2r blocks and the (p, g)-th 
block is given by D\, x , £+(0). We have: 

(43) 

^, (^p :E p+l) + (/iu ± (p)( :E p) /iu±(p+l)( a; p+l))^ 7 a;J,/u'±(p)( :r p) ( x p-l x p) + (/tu±(p-l)( :r p-l) /w± (p) (^pD^jJ, fw± (p) (^p) 

P ^/ ( x 'p— x 'p+l) 2 +(.fw±(p)( x 'p) — 7w±(p+l) ( x 'p+l)) 2 \J ( x 'p~ x 'p-l) 2 + (fw± (p) { x 'p) — fw ± (p-1) (^p-l)) 2 

A simple calculation (using (43) ) shows that all the blocks D*, x , £+(0) are zero except the 
ones with p = q,p = q + 1 and 5 = p + 1. From (43) we obtain 

fw+( P )(0)) = t-A. 



^,£ + (0)=2(I/ + W+ (p)HessA 

^..A(°) = ir 



ui+(p), 



"p+1 



□ 



In the elliptic case, det i? 2r is a polynomial in cos<x,/2 (in coshaj/2 in the hyperbolic 
case) of degree 2r(n — 1). Here, in the elliptic case, {e ±iai , ...e^"- 1 } are the eigenvalues of 
the Poincare map P 7 . 

Proposition 7.2. We have 

n-l 

detH 2r = L 2ril - n ^ Y[(2 - 2 cos raj). 

We will use Proposition 7.3 in the following subsection to prove Proposition 7.2. 

7.1.1. Poincare map and Hessian of the length functional. The linear Poincare map P 7 of 
7 is the derivative at 7(0) of the first return map to a transversal to $* at 7(0). By a non- 
degenerate periodic reflecting ray 7, we mean one whose linear Poincare map P 7 has no 
eigenvalue equal to one. For the definitions and background, we refer to [PS, KT]. 

There is an important relation between the spectrum of the Poincare map P 7 of a periodic 
m-link reflecting ray and the Hessian H m of the length functional at the corresponding 
critical point of L : (dVt) m — > R. For the following, see [KT] (Theorem 3). 

Proposition 7.3. Let 7 be a periodic m-link reflecting ray in plane domain fl. Then we 
have: 

m 

det(J - P 7 ) = det(// m )(J]// PiP+1 )- 1 , 

P =i 
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where H PiP+ i is the (p,p+ l)-th entry of H m . 

Proposition 7.3 is stated only for the plane domains. One can probably prove it for higher 
dimensions, but the formulae for the plane domains is enough for us to prove Prop. 7. 2 

Proof of Prop. 7. 2. Let us first assume n — 2. Let A r , A" 1 be the eigenvalues of P 7 r, so 
that det(J — P 7 r) = 2 — (A r + A^ 1 ). Since in our case H PjP+1 = — 1/L, from Prop. 7. 3 it follows 
that 

(44) det(J -P r ) = L 2r det H 2r ; (7 2 - link.) 



This is because if the eigenvalues of P 7 are {e } (say, in the elliptic case) then those of 
P 7 r are {e ±ira }, hence the left side of (44) equals 2 — 2cosra. Now assume n > 2 and assume 



{e 



...e 



- 1 } are the eigenvalues of P 7 . We just showed that if we define 



Hj,2r - ~T { 



l 3,+ 



1 aj- 1 ... 



1 



l 3,+ 



1 







*3r 



2rx2r 



L 2r (2 — 2 cos(raj)). Now we notice because all the blocks of the matrix 
H 2r are diagonal matrices, therefore they commute and we can write 



then det-Pj j2 r 



n-l 



det# 2 r = det(Diag(deti/i i2r , • • • , det// n _i, 2r )) = L 2r{l - n) JJ(2 - 2 cosra,-). 



We now consider the inverse Hessian H + = H^ 1 , which will be important in the calculation 
of wave invariants. We denote its matrix elements by W+ vq which corresponds to the (i, j)-th 
entry of the (p, q)-th block of the matrix % + . We also denote by %- the matrix in which 
the roles of A + ,A_ are interchanged; it is the Hessian of We also notice since B.2 r is a 
block matrix with each block a diagonal matrix so is its inverse H + . Hence the only non-zero 
entries of the inverse Hessian H + are of the form W+ vq . 

Proposition 7.4. The diagonal matrix elements h l + pp are constant when the parity of p is 
fixed, and for every 1 < i < n — 1 we have: 



p odd 



h 



ii,U 



p even 



h 



± 



h 



:,22 



± 



hi' 22 = /i!!' 11 



Proof. It is enough to show this for n — 2. This is because H 2r is a block matrix with 
commuting blocks and each block is diagonal. In fact based on our definition in (7.1.1) we 
have 
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(45) 



Hence it is enough to prove the proposition for the entries of the inverse of the 2r x 2r matrix 
H it 2r- This reduces the problem to n — 1 = 1. 

Now let us introduce the cyclic shift operator on IR 2r given by Pej = ej + i, where {e-,} 
is the standard basis, and where Pe2 r = e\. It is then easy to check that PU+P' 1 = 
hence that PU^p- 1 = HZ 1 . Since P is unitary, this says 



h p J = (H-e p ,e q ) = (PH^P-^e,) = (Hz'P'^ P~ x e q ) = h^' 1 . 

It follows that the matrix H± is invariant under even powers of the shift operator, which 
shifts the indices p — > p + 2k (k = 1, . . . , r). Hence, diagonal matrix elements of like parity 
are equal. 

□ 

7.1.2. Inverse Hessian at (Z/2Z) n - symmetric bouncing ball orbits. We first observe that in 
the case of (Z/2Z) n symmetric domains, the (2r)(n — 1) x (2r)(n— 1) Hessian of Proposition 
7.1 simplifies to: 



(46) 



2r 



-1 

T 



' A 


/ 







I 


I 


A 


I 









< 


I 


A 


/ 




> 








I 


A 


I ... 















which is a 2r x 2r block matrix, each block of size (n — 1) x {n — 1). Here A = Diag(aj) 

and / is the rank n — 1 identity matrix. The diagonal entries aj are given by 

(47) 

a, = 2cosa J /2 = — 2(l+Lvj) (elliptic case), a, = 2 cosh ctj/2 = — 2(1+Lvj) (hyperbolic case). 

We can express the inverse Hessian matrix elements h 2 J f 9 in terms of Chebychev polyno- 
mials T m , resp. U m , of the first, resp. second, kind. The Chebychev polynomials are defined 
by: 

Sill ( 771 -\- Xju 

T m (cos#) = cos m9, U m (cos9) = : — . 

Proposition 7.5. [Z3] Suppose that H 2r is given by (46). Then the matrix elements of 
are given by 

& = 2[ i-r 2r U/2)] [U2r- q+P -i(-a t /2) + U^i-ai/2)}, (1 < p < q < 2r; 1 < i < n - 1) 



h%r = o, j. 
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Proof. This formula was proved in [Z3] for the case n = 2. For the general case n > 2, we 
just use the equation (45) and reduce it to n — 2. □ 

We note that h' l ^ pq = h l ^ qp so this formula determines all of the matrix elements. It follows 
in the elliptic case that 

{-l)P-i f sm(2r-q+p)a t /2 + sin(g-p)q t /2 ^ (\ < p < q < 2r) 



(48) ti 2 l ? q = 



2(1 — cosrai)\ sinai/2 sinai/2 

(-I)'-* / sin(2r- P +g) at /2 sin( P - g )a/2 x Q < „ < „ < 2r ) 
,;) V sma,72 sin a,- / 2 / V — ^ — ^ — / 



2(1— cosrcti)\ sincti/2 sinai/2 

There are obvious analogues in the hyperbolic and mixed cases. 
The case of interest to us is 

(49) = = cot ^. 

v ; 2r 2(1 -cos r«i) sinai/2 sinai/2 2 

7.2. Properties of the phase function C + and the amplitude o^ r ' + . Since C + and 
functions of 2r(n — 1) variables £ p where 1 < p < 2r and 1 < j < n — 1, to simplify 
our notations we denote: 

Let [7] = (7^), 1 < p < 2r, 1 < j < n — 1 be a 2r x (n — 1) matrix of indices. We let 
m — I W I — 2~2 P ,j 7p- Then we define 

d m 

jym _ 

[7] ~ {dx\)^l...{dx n 2 ; 1 )^r 1 ' 

We will use 7 P for the p-th row of [7], and sometimes if 7 P , 7 g , ... are the only non-zero rows 
of [7] we write D™ p for to emphasize that 7,, 7 9 , ... are the only non-zero rows of [7]. 
The calculation of the highest derivative terms of the Balian-Bloch wave invariants uses only 
the following properties of the phase and principal amplitude which may be derived directly 
from the formulae in Theorem 6.3. 

The following Lemma is the higher dimensional generalization of Lemma 4.5 of [Z3]. It is 
proved in the same way, and the proof is therefore omitted. 

Lemma 7.6. The phase and principal amplitude of the principal oscillatory integrals 1%^ 
have the following properties: 

1. In its dependence on the boundary defining functions f±, the amplitude a 2 l' + has the 
form ai r {k,x',f±,f' ± ). 

2. As above, in its dependence on x' 

a 2r x li ■ ■ ■ 1 X 2r) = £+( X li ■ ■ ■ j x 2r)^2r (^) X l ) • ■ • J X 2r ) + ^ ~gfo^2r (^5 X l j ■ ■ • j X 2r) ) 

where 

A 2r ' (k, x 1: . . . , x 2r ) = H p=1 A p (x p , x p+1 ), 2r + l = l. (see (39)) 

3. At the critical point, the principal amplitude has the asymptotics 

<•+(£;, 0) ~ (2rL)L-^ r A(r) + 0(Ar 1 ), 
where A(r) depends only on r and not on Q. 



26 



HAMID HEZARI AND STEVE ZELDITCH 



4 - — ~ < 2ri > • 4 < r > v °^)d + o(t-» (c/.4). 



5. Vafp (fc, xi, . . . , x' 2r )\ x , =0 = 

6. D^ 1 £ + \ x , =0 = 2w + (p)D^ 1 f w+(p) (0) mod R 2r (J*f + (0), J*/_(0)). 

7. J D| +2 £ + U =0 = 2«; + (p) J D| +2 ^ + ( P )(0) mod i? 2r (^7 + (0), J*/_(0)). 

8. // [7] /ias more £/ian one non-zero row, say then 

D%£ r £ + (0) = mod R 2r (J 2 if + (0),J 2 if-(0)), 

and 

D%* r £ + {0) = mod R 2r (J^f + (0),J^f40)). 
Above, = means equality modulo lower order derivatives of f. 

7.3. Stationary phase diagrammatics. We briefly review the stationary phase expan- 
sion from the diagrammatic point of view. For more details we refer to [A, E, Z3]. 

The stationary phase expansion gives an asymptotic expansion for an oscillatory integral 



Z(k) = [ a(x)e lkS(x) dx 



where a G C£°(M. n ) and where S has a unique non-degenerate critical point in supp(a) at 
x — 0. Let us write H for the Hessian of S at 0. The stationary phase expansion is: 

Z(k) = A^ e - 7 ==e^Z A (k), 
k y/\detH\ 

where 

z A (k) = ±>r> E 'M 

3=0 (r,e)eOv,i,i-v=j y ' 

where Qvj consists of labelled graphs (T,£) with V closed vertices of valency > 3 (each 
corresponding to the phase), with one open vertex (corresponding to the amplitude), and 
with I edges. The function £ 'labels' each end of each edge of T with an index j e {1, . . . , n}. 

Above, S(T) denotes the order of the automorphism group of I, and Ie(T) denotes the 
'Feynman amplitude' associated to (r, £). By definition, Ie(r) is obtained by the following 
rule: To each edge with end labels j, k one assigns a factor of %h? k where H" 1 = (h jk ). To 
each closed vertex one assigns a factor of i dxil „; dx i v where v is the valency of the vertex and 
%x . . . , % v at the index labels of the edge ends incident on the vertex. To the open vertex, one 
assigns the factor a f 1 , where v is its valence. Then Ie(T) is the product of all these 
factors. To the empty graph one assigns the amplitude 1. In summing over (T,£) with a 
fixed graph T, one sums the product of all the factors as the indices run over {1, . . . , n}. 

We note that the power of A; in a given term with V vertices and I edges equals k~ Xr ' , 
where xv = V — I equals the Euler characteristic of the graph T' defined to be T minus 
the open vertex. We note that there are only finitely many graphs for each x because the 
valency condition forces / > 3/2 V. Thus, V < 2j, I < 3j. 
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7.4. The stationary phase calculations of I^: The data D^ 2 f±{0). In this section 
we will repeatedly use different parts of Lemma 7.6 without quoting them. 

We first claim that in the stationary phase expansion of I 2 r + P , the data D 2 i +2 /±(0) appears 
first in the k~i term . This is because any labelled graph (T,£) for which Ie(T) contains the 
factor D^i +2 /±(0) must have a closed vertex of valency > 2j + 2, or the open vertex must 
have valency > 2j + 1. The minimal absolute Euler characteristic |x(r')| in the first case is 
j. Since the Euler characteristic is calculated after the open vertex is removed, the minimal 
absolute Euler characteristic in the second case is j + 1 (there must be at least j + 1 edges.) 
Hence the latter graphs do not have minimal absolute Euler characteristic. More precisely, 
we have: 

Proposition 7.7. In the stationary phase expansion of 1®^, the only labelled graph (T,£) 
with x(r') = V — I = —j and /^(r) containing D 2 i +2 f±(0) is given by: 

• Tij+i G (i-e. V = 1,1 = j + I). The graph ^ij+i has no open vertex, one 
closed vertex and j + 1 loops at the closed vertex. 

• The only labels producing the desired data are those i v , with 1 < p < 2r fixed, which 
labels all endpoints of all edges as (i,p) where 1 < % < n — 1. (Notice the label (i,p) 
corresponds to the variable x l p .) 

In addition, the sum of the Feynman amplitudes corresponding to the labelled graphs (Tij+i, £ p ) 
above, for a fixed p, is 

$> P (r w ) = (ArL)L-^A(r)^ 2 £ MiH(fc^)* w + (p)D 2 2 £ 2 f w+(p) (0) 

l>l=i+i lp ' 

where we neglect terms with < 2j + 1 derivatives. 

Proof. We argued diagrammatically that the power /c - - 7 is the greatest power of k in which 
/±(0) appears. We also showed that a labeled graph with Euler characteristic — j which 
produces Dli +2 f±(0) must have a closed vertex of valency > 2j + 2. Now it is clear that 
such graph must have only one closed vertex and j + 1 loops. This proves the first part of 
the proposition. The second part follows easily from Lemma 7.6. 

Now let us determine J2e h p {^i,j+i) for the labelled graphs (Ti J+ i, £ p ) above. We have 
(50) 



7i+-+7? 1= j'+i 



So by Lemma 7.6 and using short-hand notations for multi-indices we get 

$> P (r 1; , +1 ) = (4rL)L-^A(ry+ 2 £ k+^( h **P)% w + (p)D 2 2 £ 2 f w+(p) (0) 



l>l=i+i 7p ' 



□ 



7.5. Proof of Theorem 2. Now we are ready to prove Theorem 2. The discussion above 
shows that modulus derivatives of order < 2j + I 

^- (2rL) L u fc) • 
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We notice that S(Tij + i) = lAutfiij+i)] = 2 j+1 (j + 1)!. We then break up the sums over 
p of even/odd parity and use Proposition 7.4 to replace the odd parity Hessian elements by 
h 1 ^ and the even ones by h 2 +. Taking into account that w + (p) = 1(— 1) if p is even (odd), 
we conclude that ( by the formula in Proposition 7.7) 

B^-JSfr £ W/+(0) - («£)'c§r/-(0)}. 

So far we have proved all parts of Theorem 2 except the last part which finds a formula for 
the wave invariants in the case of symmetries. 

7.5.1. Balian-Bloch invariants at bouncing ball orbits of (Z/2Z) n symmetric domains. Now 
if we assume the (Z/2Z) n symmetry assumptions, namely f + = f = — /_ and / being even 
in all variables, then using (49) the formula above simplifies to 

(51) 

+{polynomial of Taylor coefficients of order < 2j}. 
This finishes the proof of Theorem 2. Q.E.D. 

7.6. Recovering the Taylor Coefficients and the Proof of Theorem 1. First of all 
we prove the following lemma 

Lemma 7.8. // {ct\, . . . are linearly independent over Q then the functions 

ra x 1 
cot — 

are linearly independent over C as functions ofreN. 
Proof. Suppose that there exist coefficients such that 



c^cot^ 7 = 0, Vr eN. 



Consider the function 



ip(z u . . . , Zn-i) := ^2 c 7 ( cot • 

7 

This function is meromorphic and periodic of period 27r in each variable Zj, so it may be 
regarded as a meromorphic function on (C/Z)" -1 . It vanishes when Zj = raj/2 modulo 2ir 
for all r = 1, 2, 3, ... . But such points are dense in the real submanifold (IR/Z) n_1 and hence 
the function vanishes identically on (C/Z) n ~ 1 . This is a contradiction since the functions 
YYjZi w ^ are independent functions and by the change of variables Wj = cot Zj the functions 
YYjZi (cot ZjY j must also be independent. □ 

Now assume C I" is a domain in the class V L defined in (3). Take a non-degenerate 
bounding ball orbit 7 of length 2L which satisfies all the properties listed in (3). We would 
like to use mathematical induction and recover the Taylor coefficients of the function / where 
/ and — / are the local defining functions of dVL near the top and bottom of the bouncing ball 
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orbit respectively. First, it is possible to recover all the otj, 1 < j < n—1, under a permutation 
[Fr]. This is because | det(J — P 7 r)\ is a spectral invariant (the 0-th wave invariant). But 
we know that | det(J — P 7 r)\ = YYjZi (2 — 2 cos(raj)). Hence YYjZi (sin 2 (ro; : ,/2)) is a spectral 
invariant for all r G N. It is easy to see that this condition determines aj under a permutation. 
We fix this permutation and we argue inductively to recover all the Taylor coefficients. Since 
/ is even in all the variables, the odd order Taylor coefficients are zero. Now assume 
are given for all |7| < j. Hence the remainder polynomial term of (51) is given. Now by the 
above lemma, since all the functions (cot are linearly independent, we can recover the 
Taylor coefficients /(0). This concludes the proof of Theorem 1. 

The analogous arguments will follow in the hyperbolic or mixed hyperbolic-elliptic cases. 
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